function X = modcrop(X, modulo)
imsize = size(X);
imsize = imsize(1 : 2);
imsize = imsize - mod(imsize, modulo);
X = X(1:imsize(1), 1:imsize(2), :);
